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Abstract 

We analyze properties of superconducting state (for both s-wave and d-wave 
pairing), appearing on the "background" of the pseudogap state, induced by 
fluctuations of "dielectric" (AFM(SDW) or CDW) short - range order in 
the model of the Fermi surface with "hot spots". We present microscopic 
derivation of Ginzburg - Landau expansion, taking into account all Feynman 
diagrams of perturbation theory over electron interaction with this short - 
range order fluctuations, leading to strong electronic scattering in the vicinity 
of "hot spots" . We determine the dependence of superconducting critical tem- 
perature on the effective width of the pseudogap and on correlation length of 
short - range order fluctuations. We also find similar dependences of the main 
characteristics of such superconductor close to transition temperature. It is 
shown particularly, that specific heat discontinuity at the transition tempera- 
ture is significantly decreased in the pseudogap region of the phase diagram. 
PACS numbers: 74.20.Fg, 74.20.De 
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I. INTRODUCTION 



Pseudogap state observed in a wide region of the phase diagram of high - temperature 
superconducting cuprates is characterized by numerous anomalies of their properties, both 
in normal and superconducting states [1,2]. Apparently, the most probable scenario of pseu- 
dogap state formation in HTSC - oxides can be based [2] on the picture of strong scattering 
of current carriers by fluctuations of short - range order of "dielectric" type (e.g. antiferro- 
magnetic (AFM(SDW)) or charge density wave (CDW)) existing in this region of the phase 
diagram. In momentum space this scattering takes place in the vicinity of characteristic 
scattering vector Q = , f ) (a - lattice constant), corresponding to doubling of lattice 
period (e.g. vector of antiferromagnetism), being a "precursor" of spectrum transformation, 
appearing after the establishment of AFM(SDW) long - range order. Correspondingly, there 
appears an essentially non Fermi - liquid like renormalization of electronic spectrum in the 
vicinity of the so called "hot spots" on the Fermi surface [2]. Recently there appeared a 
number of experiments giving rather convincing evidence for precisely this scenario of pseu- 
dogap formation [3-5]. Within this picture it is possible to formulate simplified "nearly 
exactly" solvable model of the pseudogap state, describing the main properties of this state 
[2], and taking into account all Feynman diagrams of perturbation theory for scattering 
by (Gaussian) fluctuations of (pseudogap) short - range order with characteristic scattering 
vectors from the area of Q, with the width of this area defined by the appropriate correlation 
length £ [6,7]. 

Up to now the majority of theoretical papers is devoted to the studies of models of 
the pseudogap state in the normal phase for T > T c . In Refs. [8-11] we have analyzed 
superconductivity in the simplified model of the pseudogap state, based on the assumption of 
existence of "hot" (flat) patches on the Fermi surface. Within this model we have constructed 
Ginzburg - Landau expansion for different types (symmetries) of Cooper pairing [8,10] and 
also studied the main properties of superconducting state for T < T c , solving the appropriate 
Gor'kov's equations [9-11]. At first stage, we have considered greatly simplified ("toy") 
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model of Gaussian fluctuations of short - range order with an infinite correlation length, 
where it is possible to obtain an exact (analytic) solution for the pseudogap state [8,9], further 
analysis of more realistic case of finite correlation lengths was performed both for "nearly 
exactly" solvable model of Ref. [10] (assuming the self - averaging nature of superconducting 
order parameter over pseudogap fluctuations) and for a simplified exactly solvable model 
[11], where we have been able to analyze the effects due to the absence of self - averaging 



The aim of the present paper is to analyze the main properties of superconducting 
state (for different types of pairing), appearing on the "background" of the pseudogap of 
"dielectric" nature in more realistic model of "hot spots" on the Fermi surface. Here we 
shall limit ourselves to the region close to superconducting temperature T c and perform an 
analysis based on microscopic derivation of Ginzburg - Landau expansion, assuming the self 
- averaging nature of superconducting order parameter and generalizing similar approach 
used earlier for "hot patches" model in Ref. [10]. 

II. "HOT SPOTS" MODEL AND PAIRING INTERACTION. 

In the model of "nearly antiferromagnetic" Fermi - liquid, which is actively used to 
describe the microscopic nature of high - temperature superconductivity [13,14], it is usually 
assumed that the effective interaction of electrons with spin fluctuations of antiferromagnetic 
(AFM(SDW)) short - range order is of the following form: 



where g - is some interaction constant, £ - correlation length of spin fluctuations, Q = 
(n/a,n/a) - vector of antiferromagnetic ordering in dielectric phase, cu s f - characteristic 
frequency of spin fluctuations. Both dynamic spin susceptibility and effective interaction 
(1) are peaked in the region of q ~ Q), which leads to the appearance of "two types" of 
quasiparticles - "hot" one, with momenta in the vicinity of "hot spots" on the Fermi surface 



[9,11,12]. 




(1) 
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(Fig. 1), and "cold" one, with momenta close to the parts of the Fermi surface, surrounding 
diagonals of the Brillouin zone [6]. This is due to the fact, that quasiparticles from the 
vicinity of the "hot spots" are strongly scattered with the momentum transfer of the order 
of Q, due to interaction with spin fluctuations (1), while for quasiparticles with momenta 
far from these "hot spots" this interaction is relatively weak. 

For high enough temperatures ixT ^> u s f, we can neglect spin dynamics [6], limiting 
ourselves to static approximation in (1). Considerable simplification, allowing to analyze 
higher - order contributions, can be achieved by substitution of (1) by model - like static 
interaction of the following form [7]: 

« = ^ i-'X-Q-n-'X-Q.)' (2) 

where W is an effective parameter with dimension of energy. In the following, as in Refs. 
[6,7], we consider parameters W as £ phenomenological (to be determined from the exper- 
iment). Anyhow, Eq. (2) is qualitatively quite similar to the static limit of (1) and almost 
indistinguishable from it in most interesting region of |q — Q| < determining scattering 
in the vicinity of "hot spots" . 

The spectrum of "bare" (free) quasiparticles can be taken as [6] : 

£ p = — 2t(cos p x a + cos p y a) — At cos p x a cos p y a — \i (3) 

where t is the transfer integral between nearest neighbors, while t' is the transfer integral 
between second nearest neighbors on the square lattice, a is the lattice constant, ji - chemical 
potential. This expression gives rather good approximation to the results of band structure 
calculations of real HTSC - system, e.g. for YBa 2 Cu30 e+ s we have t = 0.25eV, t' = — 0.45t 
[6]. Chemical potential \x is fixed by concentration of carriers. 

The least justified is an assumption of the static nature of fluctuations, which can be 
valid only for rather high temperatures [6,7]. For low temperatures, particularly in su- 
perconducting phase, spin dynamics can become quite important, e.g. for microscopies of 
Cooper pairing in the model of "nearly antiferromagnetic" Fermi liquid [13,14]. However, 



we assume here that the static approximation is sufficient for qualitative understanding of 
the influence of pseudogap formation upon superconductivity, which will be modelled within 
phenomenological BCS - like approach. 

In the limit of infinite correlation length £ — > oo this model acquires an exact solution [15]. 
For finite £ we can construct "nearly exact" solution [7] , generalizing the one - dimensional 
approach, proposed in Ref. [16]. Then we can (approximately) sum the whole diagrammatic 
series for the one - particle electronic Green's function. 

For the contribution of an arbitrary diagram for electronic self - energy, in the iV-th 
order over the interaction (2), we write down the following Ansatz [7,16]: 

2JV-1 

^ N \e nV ) = W 2N ]J G ok] (s nP ), 

Gokj (e n p) _ ^ ^ + k ( 4 ) 

where k = kj - is the number of interaction lines, surrounding the j-th (from the 
beginning) electronic line in a given diagram, e n = 2nT{n + 1/2) (assuming e n > 0). 



Up) 



£p+q for odd k 

(5) 

£ p for even k 



Vk 



\v x (p + Q)| + \v y (p + Q)| for odd A; 

(6) 

|i>x(p)| + ky(p)| for even A; 



where v(p) = ^ - velocity of a "bare" (free) quasiparticle. 

In this approximation the contribution of an arbitrary diagram is determined, in fact, by 
the set of integers kj. Any diagram with intersection of interaction lines is actually equal to 
some diagram of the same order without intersections and the contribution of all diagrams 
with intersections can be accounted with the help of combinatorial factors s(kj) attributed 
to interaction lines on diagrams without intersections [16,7,6]. 

Combinatorial factor: 

s(k) = k (7) 
5 



in the case of commensurate fluctuations with Q = (ir/a,Tr/a) [16], if we neglect their spin 
structure [6]) (i.e. limiting ourselves to CDW - type fluctuations). Taking into account the 
spin structure of interaction within the model of "nearly antiferromagnetic" Fermi - liquid 
(spin - fermion model of Ref. [6]), we obtain more complicated combinatorics of diagrams. In 
particular, spin - conserving scattering leads to the formally commensurate combinatorics, 
while the spin - flip scattering is described by diagrams of incommensurate case ( "charged" 
random field in terms of Ref. [6]). As a result the Ansatz (4) for one - particle Green's 
function is conserved, but the combinatorial factor s(k) takes the form [6]: 



However, in the case of two - particle processes and in the analysis of the pseudogap influence 
on superconductivity, the use of spin - fermion model leads to significant complications. In 
particular, in this model there is a significant difference between vertexes, corresponding to 
spin dependent and charge interactions, and combinatorics of diagrams for spin dependent 
vertex is different from that of (8) [6]. 

In spin - fermion model [6] the spin dependent part of interaction is usually described by 
isotropic Heisenberg model. If for this interaction we assume the Ising like form, only spin - 
conserving processes remain and commensurate combinatorics of diagrams (7) remains valid 
both for one - particle Green's function and for spin and charge vertexes. For this reason, 
in the present work we shall limit ourselves to the analysis of commensurate (7) "Ising 
like" spin - fluctuations only (AFM, SDW), as well as commensurate charge fluctuations 
(CDW). Details concerning incommensurate fluctuations of CDW - type can be found in 
Refs. [7,15,16]. 

As a result we obtain the following recurrence procedure for the one - particle Green's 
function G(e n p) (continuous fraction representation) [16,7,6]: 
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(9) 



i£ n - &(p) + ikv k K - W 2 s(k + l)G k+1 (e n p) 



and "physical" Green's function is determined as G(e n p) = G (e n p). 
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Ansatz (4) for the contribution of an arbitrary diagram of iV-th order is not exact in 
general case [7]. However, for two - dimensional system we were able to show that for 
certain topologies of the Fermi surface (4) is actually an exact representation [7], in other 
cases it can be shown [7], that this approximation overestimates (in some sense) the effects 
of the finite values of correlation length £ for the given order of perturbation theory. For 
one - dimensional case, when this problem is most serious [7], the values of the density of 
states found with the help of (4) in the case of incommensurate fluctuations are in almost 
ideal quantitative correspondence [17] with the results of an exact numerical modelling of 
this problem in Refs. [18,19]. In the limit of £ — > oo Ansatz (4) reduces to an exact solution 
of Ref. [15], while in the limit of £ — > and fixed value of W it reduces to the trivial case of 
free electrons. 

This model describes non Fermi - liquid like spectral density in the vicinity of "hot 
spots" on the Fermi surface and "smooth" pseudogap in the density of states [6,7]. 

To analyze superconductivity in such a system with well developed fluctuations of short 

- range order we shall assume that superconducting pairing is determined by attractive BCS 

- like interaction (between electrons with opposite spins) of the following simplest possible 
form: 

K c (p,p / ) = -^e(p)e(p / ), (10) 

where e(p) is given by: 

1 ( s-wave pairing) 

cos{p x a) — cos(p y a) ( d x 2_ y 2-w&ve pairing) 
sin(p x a)sin(p y a) ( d xy -w&ve pairing) 
cos(p x a) + cos(p y a) ( anisotropic s-wave pairing) 
Interaction constant V is assumed, as usual, to be non zero in some interval of the width of 
2u c around the Fermi level (u> c - characteristic frequency of quanta responsible for attractive 
interaction of electrons). Then in general case the superconducting gap is anisotropic and 
given by: A(p) = Ae(p). 



e(p) = < 



In Refs. [8,9] we have analyzed the peculiarities of superconducting state in an exactly 
solvable model of the pseudogap state, induced by short - range order fluctuations with 
infinite correlation length (£ — > oo). In particular, in Ref. [9] it was shown that these 
fluctuations may lead to strong fluctuations of superconducting order parameter (energy 
gap A), which break the standard assumption of the self - averaging nature of the gap 
[20-22], which allows to perform independent averaging (over the random configurations of 
static fluctuations of short - range order) of superconducting order parameter A and different 
combinations of electronic Green's functions, entering the main equations. Usual argument 
for the possibility of such an independent averaging goes as follows [20,22]: the value of A 
significantly changes on the length scale of the order of coherence length £ ~ v F /A of 
BCS theory, while the Green's functions oscillate on much shorter length scales of the order 
of interatomic spacing. Naturally, this assumption becomes invalid with the appearance 
(in electronic system) of a new characteristic length scale £ — > oo. However, when this 
correlation length of short - range order £ <C £ (i.e. when fluctuations are correlated on the 
length scale shorter than the characteristic size of Cooper pairs), the assumption of self - 
averaging of A is apparently still valid, being broken only in case of £ > £ 1 - Thus, below we 
perform all the analysis under the usual assumption of self - averaging energy gap over the 
fluctuations of short - range order, which allows us to use standard methods of the theory 
of disordered superconductors (mean field approximation in terms of Ref. [9]). 



III. COOPER INSTABILITY. RECURRENCE PROCEDURE FOR THE VERTEX 

PART. 

It is well known that critical temperature of superconducting transition can be deter- 
mined from the equation for Cooper instability of normal phase: 



1 The absence of self - averaging of superconducting gap even for the case of £ < £oj obtained in 
Ref. [11], is apparently due to rather special model of short - range order used in this paper. 
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1-V X (0;T) = (12) 

where the generalized Cooper susceptibility is determined by diagram shown in Fig. 2 and 
equal to: 

X (q;T) = -T^^ e (p)e(p / )<f p ,p'(^,-^,q) (13) 

£n p,p' 

where $ P;P '(e n , — e n , q) is two - particle Green's function in Cooper channel, taking into 
account scattering by fluctuations of short - range order. 

Let us consider first the case of charge (CDW) fluctuations, when electron - fluctuation 
interaction is spin independent. In case of isotropic s and d XJ ,-wave pairing superconducting 
gap does not change after scattering by Q, i.e. e(p + Q) = e(p) and e(p') ~ e(p). In case 
of anisotropic s and d x 2^ y 2 pairings superconducting gap changes sign after scattering by Q, 
i.e. e(p + Q) = — e(p), thus e(p') w e(p) for p' rs p and e(p') w — e(p) for p' rs p + Q. 
Thus, for diagrams with even number of interaction lines, connecting the upper (e n ) and 
lower (—e n ) electronic lines in Fig. 2, we have p' p and the appropriate contribution to 
susceptibility is the same for both isotropic s and d xy -wave pairing. However, for diagrams 
with odd number of such interaction lines we obtain contributions differing by sign. This 
sign change can be accounted by changing the sign of interaction line, connecting the upper 
and lower lines in the loop in Fig. 2. Then for the generalized susceptibility we obtain: 

X (q; T) = -TJ2J2 G ( £ nP + q)G(-e n , -p)e 2 (p)r ± (5„, -e n , q) (14) 

Sn P 

where T ± (e n , —e n , q) is "triangular" vertex, describing electron interaction with fluctuations 
of short - range order, while the superscript ± denotes the abovementioned change of signs 
for interaction lines connecting the upper and lower electronic lines. 

Consider now the case of scattering by spin (AFM(SDW)) fluctuations. In this case 
interaction line, describing the longitudinal S z part of interaction, surrounding the vertex 
changing spin direction, should be attributed an extra factor of (—1) [6]. Because of this 
factor, in case of electron interaction with spin fluctuations contributions to generalized 
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susceptibility for different types of pairing considered above just "change places" 2 and 
susceptibility for the case of isotropic s and o^-wave pairing is determined by "triangular" 
vertex T~ , while for anisotropic s and d^^-wave case by "triangular" vertex T + . 

Thus we come to the problem of calculation of "triangular" vertex parts, describing inter- 
action with "dielectric" (pseudogap) fluctuations. For one - dimensional case and a similar 
problem (for real frequencies, T = 0) the appropriate recursion procedure was formulated, 
for the first time, in Ref. [23]. For a two - dimensional model of the pseudogap with "hot 
spots" on the Fermi surface the generalization of this procedure was performed in Ref. [24] 
and used to calculate optical conductivity. All the details of appropriate derivation can also 
be found there. The generalization to the case of Matsubara frequencies, of interest to us 
here, can be made directly. Below, for definiteness we assume e n > 0. Finally, for "trian- 
gular" vertex we obtain the recurrence procedure, described by diagrams of Fig. 3 (where 
the wavy line denotes interaction with pseudogap fluctuations), and having the following 
analytic form: 

Ffc_ 1 (e n , — e n , q) = 

= i ± w^)G A {i + 2tc ^ kq _J^ 1)(Gk+i _ Gk J r^, -,.,<> (!5) 

where G k = G k (e n p + q) and G k = G k (— e n , — p) are calculated according to (9), v k is 
defined by (6), and are given by: 

v(p + Q) for odd k 

(16) 

v(p) for even k 



"Physical" vertex is defined as r ± (e ri , — e n , q) = (e n , —e n , q). 

To determine T c we need vertices at q = 0. Then G k = G* k and vertex parts and 
become real significantly simplifying our procedures (15). For ImG k and ReG k we obtain 
the following system of recurrence equations: 



2 This is due to the fact of spin projections change in the vertex, describing "interaction" with 
superconducting gap (restricting to the case of singlet pairing). 
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ImG _ g« + kv k n - W 2 s(k + Y)ImG k+ i 

k j-y 

ReGk = _M±gWg»i! (17) 

where D fe = (f fc (p) + W 2 s(k + l^eG^+i) 2 + (e n + A^/t - Vk 2 s(£; + 1) JmGfc+i) 2 , and vertex 
parts at q = are determined by: 

r "='^' s W i^|^ r; (18) 

Going to numerical calculations we have to define characteristic energy (temperature) 
scale associated with superconductivity in the absence of pseudogap fluctuations (W = 0). 
In this case the equation for superconducting critical temperature T c0 becomes standard 
BCS - like (for the general case of anisotropic pairing): 

1 = ^lf dPi f dPt J^L (19) 

71 n^O J ° J ° " + E « 



where rh = [ 2 ^y ] is the dimensionless cutoff for the sum over Matsubara frequencies. All 
calculations were performed for typical quasiparticle spectrum in HTSC given by (3) with 
[i = —l.3t and t'/t = —0.4. Choosing, rather arbitrarily, the value of u c = OAt and 
T c0 = O.Olt we can easily determine the appropriate values of pairing constant V in (19), 
giving this value of T cQ for different types of pairing listed in (11). In particular, for the usual 
isotropic s-wave pairing we obtain = 1, while for d x 2_ y 2 pairing we get ^ = 0.55. For 
other types of pairing from (11) the values of pairing constant for this choice of parameters 
are found to be unrealistically large and we shall not quote the the results of numerical 
calculations for these symmetries 3 . 



3 Of course, our description, based on BCS equations of weak coupling theory, is not pretending 
to be realistic also for the cases of s-wave and d x 2_ y 2 -wave pairing. We only need to define 
characteristic scale of T c q and express all temperatures below in units of this temperature, assuming 
certain universality of all dependences with respect to this scale. 
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Typical results of numerical calculations of superconducting transition temperature T c 
for the system with pseudogap obtained using our recursion relations directly from (12) are 
shown in Figs. 4,5. We can see that in all cases pseudogap ("dielectric") fluctuations lead 
to significant reduction of superconducting transition temperature. This reduction stronger 
in case of d^^-wave pairing than for isotropic s-wave case. At the same time, reduction 
of the correlation length £ (growth of k) of pseudogap fluctuations leads to the growth of 
T c . These results are similar to those obtained earlier in the "hot patches" model [8,10]. 
However, significant qualitative differences also appear. From Fig. 4 it is seen, that for 
the case of isotropic s-wave pairing and scattering by charge (CDW) fluctuations, as well 
as for d x 2_ y 2 pairing and scattering by spin (AFM(SDW)) fluctuations 4 (i.e. in cases when 
the upper sign is "operational" in Eqs. (15) and (18), leading to recursion procedure for 
the vertex with same signs) there appears characteristic plateau in the dependence of T c 
on the width of the pseudogap W in the region of W < 10T c0 , while significant suppres- 
sion of T c takes place on the scale of W ~ 50T c0 . Qualitative differences appear also for 
the case of s-wave pairing and scattering by spin (AFM(CDW)) fluctuations, as well as for 
the case of d x 2_ y 2 pairing and scattering by charge fluctuations. From Fig. 5 we can see 
that in this case (when lower sign is "operational" in Eqs. (15) and (18), i.e. when we 
have recursion procedure for the vertex with alternating signs) reduction of T c is an order 
of magnitude faster. In the case of d x 2_ y2 pairing, for the values of W/T c0 corresponding 
to almost complete suppression of superconductivity, our numerical procedures become un- 
reliable. In particular, we can observe here characteristic non single valued dependence of 
T c on W, which may signify the existence of a narrow region of phase diagram with "reen- 
trant" superconductivity 5 . This behavior of T c resembles similar dependences, appearing in 



4 This last case is apparently realized in copper oxides. 

5 Our calculations show that manifestations of such behavior of T c become stronger for the case 
of scattering by incommensurate pseudogap fluctuations. 
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superconductors with Kondo impurities [25]. Alternative possibility is the appearance here 
of the region, where superconducting transition becomes first order, similarly to situation in 
superconductors with strong paramagnetic effect in external magnetic field [26]. However, 
it should be stressed that our calculations mostly show the probable appearance of some 
critical value of W/T c0 , corresponding to the complete suppression of superconductivity. In 
any case, the observed behavior requires special studies in the future and below we only give 
the results, corresponding to the region of single - valued dependences of T c . 

IV. GINZBURG LANDAU EXPANSION. 

In Ref. [8] Ginzburg - Landau expansion was derived for an exactly solvable model of 
the pseudogap with infinite correlation length of short - range order fluctuations. In Ref. 
[10] these results were extended for the case of finite correlation lengths. In these papers the 
analysis was, in fact, done only for the case of charge fluctuations and simplified model of 
pseudogap state, based on the picture of "hot" (flat) parts (patches) on the Fermi surface. 
Also in this model the signs of superconducting gap after the transfer by vector Q were 
assumed to be the same, both for s-wave and o?-wave pairings [10]. Here we shall make 
appropriate generalization for the present more realistic model of "hot spots" on the Fermi 
surface. 

Ginzburg - Landau expansion for the difference of free energy density of superconducting 
and normal states can be written in the usual form: 

F s - F n = A|A q | 2 + g 2 C|A q | 2 + ||A q | 4 , (20) 

where A q is the amplitude of the Fourier component of the order parameter, which can be 
written for different types of pairing as: A(p,q) = A ? e(p). In fact (refGiLa) is determined 
by diagrams of loop expansion of free energy in the field of random fluctuation of the order 
parameter (denoted by dashed lines) with small wave vector q [8], shown in Fig. 6. 
Coefficients of Ginzburg - Landau expansion can be conveniently expressed as: 
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A = A K A ; C = C K C ; B = B K B , (21) 

where A , Co and B denote expressions of these coefficients (derived in the Appendix) in 
the absence of pseudogap fluctuations (W = 0) for the case of an arbitrary quasiparticle 
spectrum £ p and different types of pairing: 

Ao = ^,(0)^^ < e 2 (p) >; C = N Q {0)^L < |v(p)|V(p) >; 

B o = No(0)^<e\p)>, (22) 

where angular brackets denote the usual averaging over the Fermi surface: < ... >= 
jj^oj J2 P $(£,p) • • •) an d N o (0) is the density of states at the Fermi level for free electrons. 

Now all peculiarities of the model under consideration, due to the appearance of the 
pseudogap, are contained within dimensionless coefficients Ka, Kq and Kb- In the absence 
of pseudogap fluctuations all these coefficients are equal to 1. 

Coefficients K A and K C) according to Fig. 6(a) are completely determined by the gen- 
eralized Cooper susceptibility [8,10] x(q;T), shown in Fig. 2: 

Ka = mn^m (23) 

Ao 

Generalized susceptibility, as was shown above, can be found from (14), where "triangular" 
vertices are determined by recurrence procedures (15), allowing to perform direct numerical 
calculations of the coefficients K A and K C - 

Situation with coefficient B is, in general case, more complicated. Significant simplifi- 
cations arise if we limit ourselves (in the order |A 9 | 4 ), as usual, to the case of q = 0, and 
define the coefficient B by diagram shown in Fig. 6(b). Then for coefficient Kb we get: 



K " = ^ E E e\p)(G(e nP )G(-e n , -p)) 2 (r ± ( £ „, -e n , 0)) 4 (25) 



It should be noted that Eq. (25) immediately leads to positively defined coefficient B. This 
is clear from G{— e n , — p) = G*(e n p), so that G(e n p)G(— e n , — p) is real and accordingly 
r ± (e n , — e n ,0), defined by recurrence procedure (18), is also real. 
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V. PHYSICAL PROPERTIES OF SUPERCONDUCTORS IN THE PSEUDOGAP 

STATE. 



It is well known that Ginzburg - Landau equations define two characteristic lengths - 
coherence length and penetration depth. Coherence length at a given temperature £(T) 
determines characteristic length scale of inhomogeneities of the order parameter A: 



In the absence of the pseudogap: 



em = -5. (26) 



&cs(T) = (27) 



Thus, in our model: 



em k, 



c 



e BCS (T) k a 

For the penetration depth of magnetic field we have: 



(28) 



A 2 CO = -J-^ (29) 



327re 2 AC 
Then, analogously to (28), we obtain: 

m ( K » V /2 (30) 

\bcs(T) \K a K c ) ' 1 ' 

Close to T c the upper critical field H c2 is defined via Ginzburg - Landau coefficients as: 

Hc2 ~2ne(T)~ 2vrC" ^ ' 

where 0o = cn/\e\ is the magnetic flux quantum. Then the derivative (slope) of the upper 
critical field close to T c is given by: 

dH c2 



dT 



167T0Q < e 2 (p) > Ka ,„ 2) 
7C(3)<|v(p)|V(p)> C K C 1 } 



Specific heat discontinuity at the transition point is defined as: 
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{Cs - = § (t=t) ' (33) 

where C s and C n are specific heats of superconducting and normal states. At temperature 
T c0 (in the absence of the pseudogap, W = 0): 

(C. - C n ) Tc0 - 7V(0) ?c(3) < g4(p) > . (34) 

Then the relative specific heat discontinuity in our model can be written as: 

^ _ (C s - C n ) Tc _ T c K\ 
(C s — C n )r c0 T c0 K B 

Numerical calculations of Ka, K b , K c were performed for the same typical parameters of 
the model as calculation of T c described above. The values of these coefficients, are not 
of great interest and we drop these data 6 presenting in Figs. 7-12 appropriate W/T cQ - 
dependences of physical characteristics, defined by Eqs. (26) - (35). In accord with typical 
behavior of T c described above, here we can also observe two qualitatively different types of 
behavior cue to either constant or alternating signs in recurrence equations for the vertex 
part (upper or lower signs in (15), charge or spin fluctuations). Results of our calculations of 
physical characteristics for the first type of behavior (s-wave pairing and scattering by charge 
(CDW) fluctuations, d x 2_ y 2-wave pairing and scattering by spin (AFM(SDW)) fluctuations) 
are shown in Figs. 7 - 10. We can see that with the growth of pseudogap width W 
coherence length £(T) is reduced, while penetration depth A(T) grows, as compared to the 
appropriate values of BCS theory. Dependence of both characteristic lengths on parameter 
na is rather weak and in Figs. 7, 8 we present data for only one value of Ka = 0.2. The slope 
(derivative) of the upper critical field at T = T c initially grows, then starts to diminish. 
Most characteristic is the suppression of specific heat discontinuity (as compared to BCS 
value) shown in Fig. 10, which is in direct qualitative accord with experimental data [27]. 



6 Typical dependences of these coefficients on the parameter W/T c q are functions rather rapidly 
decreasing from 1 in the region of existence of superconducting state. 
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Note that for specific heat discontinuity we also have noticeable plateau in the region of 
W/T c0 < 10, analogous to that observed above in similar dependence of T c . 

Behavior of physical characteristics of superconductor for the case of s-wave pairing and 
scattering by spin (AFM(CDW)) fluctuations as well as for d x 2_ y 2 pairing and scattering by 
charge (CDW) fluctuations are shown in Figs. 11, 12. Data on characteristic lengths are 
not given, as both coherence length £(T), and penetration depth A(T) are practically the 
same as in BCS theory in almost all region of existence of superconducting state (except 
very narrow region where T c can be non single valued and disappear, where both lengths are 
rapidly growing). As to the slope of an upper critical field and specific heat discontinuity 
at superconducting transition, both are rapidly suppressed with the growth of W/T c0 , up 
to the region of non single valued T c (the region of possible "reentrant" behavior or type I 
transition) (cf. Fig. 11, 12), where both are characterized by non monotonous behavior. 

VI. CONCLUSION 

In this paper we have analyzed peculiarities of superconducting state forming on the 
"background" of the pseudogap induced by electron scattering by fluctuations of "dielectric" 
short - range order in the model of "hot spots" on the Fermi surface. Our analysis is based on 
microscopic derivation of Ginzburg - Landau expansion and takes into account all high order 
contributions of perturbation theory for scattering by pseudogap fluctuations. "Condensed" 
phase of such superconductor should be described by appropriate Gor'kov's equations for a 
superconductor with the pseudogap [10], which will be the task for the future work. 

Our main conclusion is that superconductivity is suppressed by pseudogap fluctuations of 
CDW, AFM(SDW) type with two qualitatively different model of such suppression, related 
to either constant or alternating signs in recurrence equations for the vertex part (upper 
or lower signs in (15, charge or spin fluctuations). While the case of scattering by spin 
fluctuations and pairing with d x i_ y i - symmetry is apparently realized in high - temper- 
ature superconducting copper oxides, we are not aware of systems with unusual behavior 
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obtained above for the case of s-wave pairing and scattering by spin (AFM(CDW)) fluctua- 
tions, as well as for the case of d x 2_ y 2 pairing and scattering by charge (CDW) fluctuations. 
Experimental search of systems with this behavior may be of significant interest. 

The most important problem in the description of the pseudogap state in real HTSC - 
cuprates is an explanation of doping dependences of the main physical characteristics. In 
our model, these dependences should be expressed via corresponding doping dependences of 
the effective pseudogap width W and correlation length £. Unfortunately, these dependences 
are determined from rather indirect experimental data and are poorly defined [1,2] 7 . As a 
very crude statement, we can claim that correlation length £ does not change significantly 
in rather wide region of doping concentrations, while the pseudogap width W drops linearly 
with concentration growth from the values of the order of 10 3 K close to the insulating 
phase, to the values of the order of T c at optimal doping, becoming zero at slightly higher 
concentrations of current carriers (Cf. Fig.6 in the review paper [2], which is based of Fig. 4 
of Ref. [3], with some compilation of appropriate data for YBCO). Using such doping 
dependence we can immediately reexpress W - dependences given above as dependences on 
concentration of doping impurity. In the oversimplfied version of our model with an infinite 
correlation length and Fermi surface with complete nesting, assuming also linear doping 
dependence for T c0 , such calculations were given in a recent paper [28]. It was shown that 
even within this simplest approach the qualitative form of the phase diagram of cuprates can 
be nicely reproduced. However, the obvious crudeness of the model, as well as the absence 
of reliable data on concentration dependences of W, £ and T c0 , do not allow any serious 
"improvement" of these qualitative conclusions. 

Among shortcomings of the present model, besides many times mentioned neglect of 
dynamics of short - range order fluctuations and limitation to Gaussian statistics of these 



Concentartion dependence of T c q may also be quite important, but in fact it is completely 
unknown. 
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fluctuations, we also note the simplified analysis of spin structure of interaction, which 
assumed its Ising like nature. It will be of great interest to perform similar analysis for the 
general case of Heisenberg interaction. 

This work is partly supported by RFBR grant 02-02-16031, as well as by the program 
of fundamental research of Presidium of the RAS "Quantum macrophysics" and program of 
the Division of Physical Sciences of the RAS "Strongly correlated electrons in semiconduc- 
tors, metals, superconductors and magnetic materials", and by the research project of the 
Russian Ministry of Science and Industries "The studies of collective and quantum effects 
in condensed matter". 
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APPENDIX A: GINZBURG LANDAU COEFFICIENTS FOR ANISOTROPIC 
SUPERCONDUCTOR IN THE ABSENCE OF THE PSEUDOGAP. 



In the absence of pseudogap fluctuations (W = 0) the generalized Cooper susceptibility 
define by diagram of Fig. 2 takes the following form: 

Xo (q;T) = -TEE^(P) , \ ■ * (Al) 

e n P n Sp+q '■^n Sp 

Then for susceptibility at q = 0, defining the coefficient A , we get: 

Xo (0;T) = -T££e 2 (p)^-- = -T W £ " ^(p) - 

- "iVo(0)TE r ^^f^ = Xbcs(0;T) < e 2 (p) > (A2) 

where angular brackets denote averaging over the Fermi surface and we introduced the 
standard susceptibility of BCS model for isotropic s-wave pairing Xbcs(®',T). 
As a result, we obtain coefficient A as: 

A = Xo (0; T) - Xo (0; T c ) = A BCS < e 2 (p) > (A3) 

where 

Abcs = Xbcs(0; T) - X bcs(0; T c ) = N (0)^—^ (A4) 

i- c 



is the standard expression for coefficient A for the case of isotropic s-wave pairing. 

Coefficient C of Ginzburg - Landau expansion is defined by susceptibility (Al) at small 

q: 



Xo(g;T c ) - Xo(0;T c ) 



Co = fen 2 *wcj (A5) 



Expanding expression (Al) for Xo(q; T c ) in powers of q we get: 

Xo(q;T c ) =Xo(0;Tj+T c EEi§fSie 2 ( P )(v( P )q) 2 (A6) 
so that coefficient C is found to be: 
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Co = T c Y, / de 4 " ^ 3 E^-e P )e 2 (p)|v( P )|W(0) 



.,,..-00 s 4(e*+£ 2 ) 3 P 
34 - e 2 



/oo \F r 

dC j " ^ 3 E^(Ue 2 (p)|v(p)|W(0) = 
-00 4(£ n -+- ? J p 



= ^o(O) < e 2 ( P )|v(p)| W(0) > (A7) 



■ c 

^00 1 



where <p is an angle between vectors v(p) and q, £(3) = Y^=i ^3 ~ 1, 202 

For square lattice both the Fermi surface and |v(p)| are symmetric with respect to 
rotation over an angle |, the same symmetry is valid for e(p) (for all types of pairing, 
mentioned above). Thus we easily find: 

< e 2 (p)|v(p)| W(0) >= \ < e 2 (p)|v(p)| 2 (l + cos(20)) >= \ < e 2 (p)|v(p)| 2 > (A8) 

because cos (2<p) changes sign after the rotation of p by an angle |. In fact, the direction of 
the velocity v(p) after this rotation changes to a perpendicular one, accordingly cos(20) — > 
—cos (20). As a result we obtain for the coefficient Co an isotropic expression: 

Co = iV (0)^^<|v(p)|V(p)>, (A9) 

which for the case of isotropic s-wave pairing and spherical Fermi surface becomes the 
standard one: 

Cbcs = ^o(0)|^ (A10) 

Coefficient B defined by the diagram shown in Fig. 6(b), in the absence of pseudogap 
fluctuations (W = 0) and for q = takes the following form: 

_ I _ poo ! 

Bo = T c £ £ 7^T7^ e4 (P) = T < E y Gl^Tw E ^ - ^(P) - 

» wg£ ^ gE gr = Bbcs < c4(p) > (A11) 



where 



= ^o(O)^ (A12) 



is the standard expression for B in case of isotropic s-wave pairing. 
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FIGURES 




FIG. 1. Fermi surface with "hot spots" connected by scattering vector of the order of 

Q = (M)- 
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FIG. 3. Recurrence equations for the vertex part. 
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FIG. 4. Dependence of superconducting transition temperature T c /T c q on the effective width 
of the pseudogap W/T c q for s-wave pairing and scattering by charge (CDW) fluctuations (curves 
si and s2) and for d^^-wave pairing and scattering by spin (AFM(SDW)) fluctuations (curves 
dl and d2). Data are given for the values of inverse correlation length na = 0.2 (si and dl) and 
kq = 0.5 (s2 and d2). 
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FIG. 5. Dependence of superconducting transition temperature T c /T c q on the effective width of 
the pseudogap W/T c q for s-wave pairing and scattering by spin (AFM(SDW)) fluctuations (curves 
si and s2) and for d x 2_ y 2-wsve pairing and scattering by charge (CDW) fluctuations (curves dl 
and d2). Data are given for the values of inverse correlation length kci = 0.2 (si and dl) and 
na = 1.0 (s2 and d2). 
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20 W/TCO 40 60 

FIG. 7. Dependence of the square of the coherence length £ 2 /£bcs on the effective width of the 
pseudogap W/T c q for s-wave pairing and scattering by charge (CDW) fluctuations (full curve) and 
for d a ,2_ J/ 2-wave pairing and scattering by spin (AFM(SDW)) fluctuations (dotted curve). Data are 
given for the value of inverse correlation length na = 0.2. 
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20 w/Tco 40 60 

FIG. 8. Dependence of penetration depth \/Xbcs on the effective width of the pseudogap 

W/T c q for s-wave pairing and scattering by charge (CDW) fluctuations (full curve) and for 
d x 2_ y 2-wave pairing and scattering by spin (AFM(SDW)) fluctuations (dashed curve). Data are 
given for the value of inverse correlation length kq = 0.2. 
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20 W/TCO 40 60 

FIG. 9. Dependence of the derivative (slope) of the upper critical field (normalized by its 
value for W = 0) on the effective width of the pseudogap for s-wave pairing and scattering by 
charge (CDW) fluctuations (curves si and s2) and for d x 2_ y 2-wave pairing and scattering on spin 
(AFM(SDW)) fluctuations (curves dl and d2). Data are given for the values of inverse correlation 
length hm = 0.2 (si and dl) and na = 0.5 (s2 and d2). 
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FIG. 10. Dependence of specific heat discontinuity at superconducting transition on the effec- 
tive width of the pseudogap W/T c q for s-wave pairing and scattering by charge (CDW) fluctuations 
(curves si and s2) and for d x 2_ y 2-wave pairing and scattering by spin (AFM(SDW)) fluctuations 
(curves dl and d2). Data are given for the values of inverse correlation length na = 0.2 (si and 
dl) and na = 0.5 (s2 and d2). 
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FIG. 11. Dependence of the derivative (slope) of the upper critical field (normalized by its value 
for W = 0) on the effective width of the pseudogap W/T c q for s-wave pairing and scattering by 
spin (AFM(CDW)) fluctuations (curves si and s2) and for d^^-wave pairing and scattering by 
charge (CDW) fluctuations (curves dl and d2). Data are given for the values of inverse correlation 
length hm = 0.2 (si and dl) and na = 1.0 (s2 and d2). 



32 



AC 

"1 1 1 1 1 1 1 1 1 1 r 




U.U1 . 1 . 1 . 1 . 1 . 1 1 
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FIG. 12. Dependence of specific heat discontinuity at superconducting transition on the ef- 
fective width of the pseudogap W/T c q for s-wave pairing and scattering by spin (AFM(CDW)) 
fluctuations (curves si and s2) and for d x 2_ y 2-w&ve pairing and scattering on charge (CDW) fluc- 
tuations (curves dl and d2). Data are given for the values of inverse correlation length na = 0.2 
(si and dl) and na = 1.0 (s2 and d2). 
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